Detecting and characterizing frequency fluctuations of vibrational modes 
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I. INTRODUCTION 

In recent years, there has been made significant 
progress in developing micro- and nano-mechanical sys- 
tems that display slowly decaying vibrations. For dif- 
ferent types of such systems, the ratio of the vibration 
eigenfrequency to the decay rate, the quality factor Q, 
has reached > lO^.'illS This has allowed studying new 
physics, including quantum phenomena,!^ and opened a 
way for numerous applications, like highly-sensitive mass 
sensing^"^ and, potentially, high-accuracy nanomechani- 
cal clocks. In parallel, high-Q modes of superconducting 
cavities have been used for control and measurement of 
Josephson-j unction based qubits.'Sl 

An important problem in the studies of nanomechan- 
ical vibrations and superconducting cavity modes is to 
understand the mechanisms of their decay and loss of 
coherence. Often one separates decay and fluctuations 
of the vibration amplitude and fluctuations of the vi- 
bration phase. Phase fluctuations are not only inter- 
esting on their own but are particularly important for 
applications, as they can impose limits on the sensitiv- 
ity of a device. They can come from the thermal noise 
that accompanies vibration decay and is a consequence 
of coupling to a thermal reservoir. A more delicate and 
often more important source is fluctuations of the vi- 
bration frequency. They can have various origins, see 
Ref. [To] and papers cited therein, with recent examples 
being random attachment or detac hment of molecules 
to a resonator that changes its m ass! '' * ^ * ^"'^ ! -'^'^ !"^ molecule 
diff'usion along the resonatorj ^"^ ! ^^ ! coupling of the vibra- 
tional mode to two-state fluctuators,^® and, for nonlinear 
vibrations, frequency modulation by fluctuations of the 
vibration amplitude.'I^ 

In this paper, we suggest a simple way of separating 
and characterizing frequency fluctuations in vibrational 
systems. In two-level systems, frequency fluctuations 
lead to the difference between the Ti and T2 relaxation 
times and are routinely separated from decay using non- 
linear response to an external field.l^ In contrast, the 
response of linear vibrations is inherently linear, and the 
spectrum of the response remains a major source of in- 
formation about the dynamics. If frequency fluctuations 
are the dominating factor, this spectrum reveals some 



of their features.'iSiniiMS! However, in many cases of in- 
terest it does not provide enough information, and often 
does not allow one to even detect frequency fluctuations 
at all. For example, for broadband Gaussian frequency 
noise, the absorption spectrum is Lorentzian, as if there 
were no frequency noise, even though the overall width 
of the spectrum exceeds the width due to decay. 

We show below that frequency fluctuations can be 
studied by using, in a different way, essentially the same 
measurement as that used to find the absorption spec- 
trum, i.e., by looking at the response of a resonantly 
modulated oscillator. This applies to both classical and 
quantum oscillators. The idea is to study such correla- 
tors of the quadrature and in-phase components of the 
oscillator displacement that are specifically sensitive to 
frequency fluctuations. As we show, these are correlators 
and moments of the complex vibration amplitude. They 
allow one not only to reveal frequency noise, but also to 
study its statistics, for both classical and quantum vibra- 
tions. The sensitivity to the noise statistics is illustrated 
for important examples of the noise. 

II. EQUATION OF MOTION FOR THE 
QUADRATURES 

We assume that the oscillator energy relaxation comes 
from the coupling to a thermal reservoir, which is linear 
in the oscillator coordinate and momentum and weak, 
so that the decay rate F <C wq, where loq is the oscil- 
lator eigenfrequency in the absence of frequency noise. 
With the noise the frequency becomes u)q +£,{!)■ We as- 
sume that the frequency noise ^(t) is a stationary process 
and {£,{!)) = 0. Of primary interest is the case of small 
noise, where overall frequency fluctuations are small com- 
pared to LOq. 

Phenomenologically, the motion of the oscillator with 
coordinate q and with unit mass in the presence of a 
driving force F cos LOpt is described by equation 

q + 2Tq+ [cjg + 2uja^{t)]q = Fcosajpt + f{t), (1) 

where f{t) is zero-mean additive thermal noise. We will 
be considering this motion in the rotating frame on times 
much longer than (jJq^ . On this scale, the approxima- 
tion of Markovian relaxation of the oscillator amplitude 
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and phase applies even where Eq. ([T]) does not apply; in 
other words, the assumption of Ohmic dissipation is not 
requiredPSl 

We assume that the typical frequencies of the noise 
S,(t) are small compared to ujq. This is the case for many 
systems of current interest; coupling of the oscillator to 
the source of such noise does not lead to energy relax- 
ation via nonlinear friction'^^ and to random parametric 
excitation of the oscillator j^^iSS, 

We consider resonant driving, with frequency lof close 
to Wq, i.e., |i5u;| <C uiq, where Soj — ujp ~ luq. The os- 
cillator dynamics can be then conveniently analyzed by 
changing from the coordinate q and momentum p — q to 
slowly varying on the time scale cop^ complex variables 
u{t),u*{t), 

q{t) = uexp(iujpt) + u* exp(—iujpt), (2) 
p(t) — iujp [uexp(iwi?t) — u* exp(— iwi?^)] . 

Function u{t) is the complex vibration amplitude. From 
Eq. Re u and Im u give, respectively, the in-phase 
and quadrature components of the oscillator displace- 
ment at frequency ujp. 

On the time scale that largely exceeds the correlation 
time of the thermal reservoir and ujp^^ both a micr oscopic 
theory for weak oscillator-to-reservoir couplinglZBSEll 
and the phenomenological model of Eq. ([!]) lead to an 
equation of motion for u of the form 



-[T + iSuj - i£,{t)]i 



iF 
Aujp 



(3) 



where fu{i) = —(i/'2'-^F)f{t)exp{—iujpt); the renormal- 
ization of luq due to the coupling to a thermal reservoir 
that emerges in the microscopic theory has been incor- 
porated into ujQ. 

For large times compared to the decay time F^^, the 
initial state of the oscillator is "forgotten" and u{t) be- 
comes a linear superposition of the terms that describe 
forced oscillator vibrations and thermal fluctuations, 



u{t) = {F/4ujp)up{t)+UtY,{t), 



(4) 



Up{t) 



I dtix* {t ~ ti) exp i I dt[^{t[ 

J —OO L ti 



= iexp(-A*i), A = A{ujp) = T + iSuj. 

Here, upit) is the scaled complex amplitude of forced vi- 
brations; x(i) is the scaled oscillator susceptibility in the 
absence of frequency fluctuations. The term uth comes 
from the additive thermal noise. 



i / diix*(i-ii)/u(ii)exp i ( dt[^{t[ 

J —OO L -^ii 



(5) 

We are interested in the effects of frequency noise, not 
additive noise. By increasing the field F, the term cx 
UF{t) in u{t), Eq. Q, can be made larger than the typical 
value of Wth- However, as we will show, the contribution 
of additive noise to the correlators of u{t) vanishes in the 



approximation used to derive Eq. (|3| , which allows one to 
use even comparatively weak driving fields for studying 
frequency noise. 



III. CORRELATORS OF THE COMPLEX 
AMPLITUDE: INDEPENDENCE OF ADDITIVE 
NOISE 

The two noises that determine the oscillator dynamics, 
/(<) and £,{t), are uncorrelated, generally. The noise f{t) 
results from the linear in q,p coupling to a thermal reser- 
voir. The major effect on the oscillator comes from the 
Fourier components of f{t) with frequencies w such that 
\uj — uJo\<^ujo- It was assumed in deriving the Markovian 
equation of motion for u{t), Eq. ([s]), that the spectral 
density of f(t) is smooth around ujq and ca n be set equal 



to a constant for 



< 



|(Su;|,(e2(t))i/2.E3Mlll 



In contrast, noise ^(i) comes from either an exter- 
nal nonequilibrium source or from the interaction with 
a thermal reservoir that is effectively quadratic in q,p 
and couples the oscillator to the degrees of freedom other 
than those that lead to f{t). 

An important consequence of the statistical indepen- 
dence of f{t) and ^{t) is that 



«h(i)> - 0, 



n = 1,2,3, 



(6) 



A simple way to see this is by comparing the expressions 
for Mth(0 and uth{t + to) with an arbitrary ^o- If one 
writes Eq. ( 5| for Mtii(i+to) and changes from integrating 
over ti and t'^ to ti — ti — to and t[ = t'^ — to, respec- 
tively, the expression for uth(t -f- to) becomes of the same 
form as Uth{t), except that the noises f{t) and ^(t) are 
evaluated for the time shifted by to and there emerges an 
extra factor exp{—iujpto) from the interrelation between 
and f{t). For stationary noises f{t) and ^(t), the 
change of the origin of time does not affect any average 
values, and therefore the only difference between averag- 
ing Mth(i) and uth(t + to) is the factor exp{—iu!pto). The 
average moments of uth{t) may not depend on to, and 
therefore they are equal to zero. So are also correlators 
(uth(ii) • . .uth(i„)). 
From the above arguments, 

(w(<i)... = iF/4u;pr{upiti)...up{t„)), (7) 

and below we will be interested in calculating correlators 
of Up. They are independent from additive noise, for a 
linear oscillator, and therefore measuring them immedi- 
ately reveals frequency noise. We note that the nth-order 
correlator is proportional to the nth power of the driving 
field; still the oscillator response remains linear. 



IV. QUASIWHITE FREQUENCY NOISE 

The explicit expressions Q and Q allow one to an- 
alyze correlators of the complex amplitude for various 
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types of frequency noise. We will discuss several noise 
models of interest for experiment and show how by mea- 
suring the correlators one can study the noise statistics. 

One of the most important is noise with a compara- 
tively broad frequency spectrum, which is flat up to a 
characteristic cutoff frequency Wcorr such that F, \6uj\ ^ 
Wcorr i^o- Such uoisc is effectively S correlated on a 
time scale long compared to uJq^ ■ It can come, for exam- 
ple, from quasielastic sca tterin g of phonons or other ex- 
citations off the oscillator j ^^ l ^^ l in which case it is approx- 
imately Gaussian, or from the discreteness of the electric 
current that modulates the oscillator, in which case it is 
close to Poissonian, or it can come from other processes 
and have a different statistics. For a (5-correlated noise, 
it is convenient to do the averaging in Eq. ([t]) using the 
noise characteristic functional, which can be written as 



exp 



exp 



i / dtk{t)Cit) 



dtfi{k{t)) 



(8) 



Function fJ.{k) is determined by the noise statistics; 
for zero- mean Gaussian noise of intensity D, with 
= 2DS{t~t'), and for zero-mean Poisson noise 
^(t) = gJ2n^i^~^n)—g'^ with pulse area 5 and pulse rate 
1/ we have, respectively, /i = /iq and n — ij,p, wher^^Il 

fiGik)^Dk^, ^ip{k) ^ V {l - e'^^ + ikg) . (9) 

For (5-correlated frequency noise, the response of the 
oscillator to the driving as a function of frequency de- 
tuning (5a; = a; — Wo is of the same functional form as 
without the noise. From Eqs. Q and ([s]), 

{u) = {F/AujF){uF{t)) = -i{F/4ujf){T + iSi:j)-\ 
T + Re nil), SCj^Suj + lmfi{l). (10) 

The noise leads to broadening of the spectrum of the 
response and, generally, to the shift of the oscillator fre- 
quency. Both are determined by the value of function fi 
for fc = 1. In particular, for quasiwhite noise, the incre- 
ment of the half- width is given by the noise intensity D, 
a well-known result, whereas for Poisson noise this incre- 
ment is J^(l — cosg), it oscillates with increasing g and 
increases with the pulse rate v. 



It follows from Eq. (10 1 that, from the oscillator spec- 
trum taken alone, one cannot tell whether there is fre- 
quency noise at all. However, the pair correlator of 
the complex amplitude makes it possible to identify the 
presence of the noise. A straightforward but somewhat 
lengthy calculation shows that 



,2M(l)-/i(2) 



2A-h^(2) 
X exp —(r + iSu})t 



(i>0).(ll) 



7^ (u)^. From 
Eq. (11), the variance of the complex amplitude u{t) is 



Because of the frequency noise. 



cx 2^(1) — /i(2), and is thus determined by the nonlinear- 
ity of the function /Lt(fc). In particular, for Gaussian and 
Poisson noises we have, respectively, 2/X(5(l) — /^g(2) — 
-2D and 2^p(l) - fip{2) = v[l - exp(ig)]2. The time 
decay of the pair correlator of Su{t) — u(t) — (u) is expo- 
nential, with exponent A + fi{l) = F -I- iSu. 

Not only does the pair correlator, Eq (111, allow one 



to reveal frequency noise where there are no conventional 
spectral signatures of it, but it also gives an insight into 
the noise statistics. More insights can be gained from the 
higher-order moments of u(t). By writing 



rO rti 

u^(0) dti / dt2 . 




At,+z(n + l-j) / dt'^atj) 



we obtain from Eq. ([8|) 



f-iF\ 

\4:UJF 



^ 3 = 1 



(12) 



FromEq. (n2|, by measuring the moments of the complex 
amplitude u(t), one can find function /i(A:) for all integer 
k and therefore, given that this function is analytical at 
least near the real-A: axis, find the whole ^(fc) and thus 
the full statistics of the ^-correlated noise £^{t). 



V. COMPARATIVELY WEAK FREQUENCY 
NOISE 

The presence of non-(5-correlated frequency noise can 
be usually directly seen in the spectrum of the oscillator 
response, if the noise is sufficiently strong. For example, 
the absorption spec trum may h ave a fine structure or 
become asymmetric .I^ISIloililillill 'pjjg situation is more 

complicated where the noise is comparatively weak, so 
that the shape of the spectrum is weakly distorted com- 
pared to the Lorentzian contour. We now show that the 
moments of the complex amplitude allow one to detect 
frequency noise and study its statistics even in this case. 

We will express the moments in terms of the correlators 
of the frequency noise. The lowest-order correlators of a 
stationary zero-mean noise are 



S2(c.) = (2^)-iy dte^-'mm), (13) 

S3(a;i,o;2) = (2^)-2 / diidi2e^("^*^+"^*^) (e(ti)e(t2)e(0)). 



Since f (ii) for different ti commute with each other, we 
have 

S2(w) S2(-w), 23(^1,012) = S3(a;2,a;i) 

= ^^{-uji - a;2,W2) = S3(tJi, -oji - 0J2). (14) 
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By expanding Eq. Q for up to third order in (,{t), we 
obtain 



\(0) 



dw 



A2 



are the same. In the absence of couphng to a thermal 
reservoir, the Hamiltonian of the osciUator in a resonant 
field F cos ojpt in the presence of weak classical frequency 
noise £_{t) is 



X A (3A^ +ujI+ UJ1UJ2 + i^l) /3, 



(15) 



where = —iF/{4:UjpA) is the complex amplitude 

that describes forced vibrations in the absence of fre- 
quency noise and 



(wi,a;2) 



^3(a;i,a;2) 



(16) 



It is clear from Eq. (15 1 that, to third order in ^(t), the 



effect of frequency noise on the spectrum of the oscilla- 
tor response can be described as renormalization of the 
decay rate F and the eigenfrequency cjq, and therefore 
from spectroscopic data it is hard to tell whether weak 
frequency noise is present at all. 

Frequency noise can be detected by measuring higher 
moments of the complex amplitude. Keeping only the 
second and third-order correlators of ^(t), we obtain for 
the variance of the complex amplitude 



duj 



r^2(w) 



—i J dwidw2S^'^^(wi, CJ2) 

X [2A^ + A(a;^ + uj^) + iu!i0j2{^^i + 1^2)] • (17) 

Here, we have used the symmetry properties of S2 and 
S3, Eq. (14); note that if the noise £^{t) has time-reversal 
symmetry, the term oc LUii02{itJi +1^2) in the integrand in 
Eq. (17) can be disregarded. 



For weak noise ^(t), the leading contribution to the 
variance of u comes from the second-order term oc S2. 
To reveal a nonzero third-order noise correlator, in ad- 
dition to the variance of u one should measure the third 
cumulant of u, 



-i I duJiduJ2^3^\uJi,uJ2) 



(^.)3 

X [A^ -I- A{ujf +UJ2+ 0J1UJ2) + iujiuj2{^i + ^^2)] ■ (18) 
One can see that, in the case of weak (5-correlated noise. 



Eqs. (TKH) — (12) agree with Eqs. (h5| — (18 1. However, 



Hq = h[ijJo + ^(<)]a^a — qF cosujpt, 



(19) 



where a — {2hujo)^^^^{ujoq + ip) is the lowering operator 
of the oscillator. 

The effect of coupling to a thermal reservoir can be 
conveniently analyzed in the rotating wave approxima- 
tion by going to the interaction representation with oper- 
ator U(t) = exp{—iujpa^ at). When typical noise frequen- 
cies are small compared to ujp, the resulting equation for 
the oscillator density matrix po, for a given realization of 
the frequency noise, has a familiar form 



Po = i[Suj - ^(i)][a^a, Po] 
+i[F'a^ +F'*a,po]. 



Tpo 



(20) 



Here, F' = (Shuja) ^/"^F, operator f describes oscillator 
decay, Fp = T{n + l){a^ ap — 2apa^ + pa^a) + Tn{aa^ p — 
2a^ pa + paa^), where n = [exp{fiujo/kBT) — is the 
Planck number; as in the classical analysis, the renormal- 
ization of the oscillator frequency due to the coupling to 
the thermal reservoir is assumed to have been incorpo- 
rated into Wo- We emphasize that the density matrix 
Po has not been averaged over the realizations of ^(i), it 
fluctuates in time. 

Equation ( 20 ) allows one to perform averaging over 
thermal fluctuations of the oscillator for given S,{t). It 
leads to a chain of equations for the moments a"(t) = 
Tr[a"po(0] of the operator a. 



dt 



a" = -71 [A* + 2^(i)]a" + inF'a 



(21) 



with A = F + iSuj, cf. Eq. If wc assume that the 
field F is turned on a diabatically at i — >■ —00, the initial 
condition to Eq. (21 ) is a" — >• for t — —00. 
The solution of Eq. (I^Tl) is 



a"(i) = [a{t)] 



a{t)^ F'u*p{t), 



(22) 



where up{t) is given by Eq. Therefore the averaging 
of the moments of operator a over realizations of £,{t), 
that gives the mean value (a"(i)), can be done in the 
same way as for a classical oscillator. Hence, the results 
for a linear classical oscillator immediately apply to a 
linear quantum oscillator. 



the results of this section are not limited to (5-correlated 



VII. MARKOV FREQUENCY NOISE 



VI. QUANTUM FORMULATION 

The above arguments can be immediately extended 
to the quantum regime, since the responses of quantum 
and classical harmonic oscillators to resonant modulation 



The analysis of the effect of frequency noise is sim- 
plified in the case where the noise is Markovian. Such 
noise can be continuous, as in the case of frequency fluc- 
tuations due to diffu sion of massive particles along a 
nanoresonator IHMiSI or discrete (i.e., takes on discrete 
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values), as in the case of random attachment or detach- 
ment of massive particles to a mechanical resonator-'il'^'^' 
or transitions between quantum states of a nonlincarly 
coupled vibrational mode of a trapped electron.^* In both 
cases, its probability distribution t) is described by 
equation p — Wp. Operator W is independent of time, 
for a stationary process. For a continuous process ^(i), 
is a differential operator with respect to ^; for exam- 
ple, for a diffusion process £,{t), the equation for p is the 
Fokker-Planck equation.'^ For a discrete process, W de- 
scribes transitions between different discrete values of ^, 
with appropriate transition rates. 

Instead of the density matrix pQ, which depends on 
a realization of the noise £,{t), for Markovian £,(t), it 
is convenient to introduce density matrix p{S.,t), which 
remains an operator with respect to the oscillator vari- 
ables, but also depends on ^ as a variable. For continuous 
C(t), we have_p(e,t) = {po{t)S{^ - where {...)^ 

means averaging over realizations of ^(t). For discrete 
$,(t), one should use the same definition, but with Kro- 
necker's delta instead of the 5 function. 

The equation for t) is an obvious extension of 

Eq. m, 



p = i[5uj-£\Wa,p]-tp+Wp 

+i[F'a^ +F'*a,pl p = p{i,t). (23) 

In this model, there is no backaction from the oscillator 
on the frequency noise source, and therefore operator W 
does not depend on the dynamical variables of the os- 
cillator. Then one can immediately write a system of 
equations for the moments A{n, ^, t) — Tr [a"/9(^, t)] of 
operator a, which will now be functions of the random 
variable ^ rather than functionals of ^ (t) , 

dtA{n, t) = -n[A* + ii]A{n, ^, t) + inF'A{n - 1, ^, t) 

+WA{n,^,t), {a"{t)) = f d^A{n,tt) (24) 



(for discrete-valued noise, the integral over d£, should be 
replaced by a sum). 

For n = 1, the stationary solution of Eq. (24) was dis- 



cussed earlier for several models of frequency noise.'i^'i^ 
As a function of detuning Suj, (a) gives the spectrum of 
the response of the oscillator to a resonant force. As men- 
tioned above, where the noise ^(t) is strong compared to 
r (but still weak compared to wq), it can significantly 
change the spectrum compared to the C = case, mak- 
ing it possible to detect the presence of the noise and find 
some of its characteristics. 

By studying the moments of the complex amplitude 
(a") one can extract much more information about the 
frequency noise than just from the spectrum. We note 
that one can think of functions A(n, ^, t) in the station- 
ary regime as "partial moments" of the oscillator for a 
given eigenfrequency ljq + Equation (24 1 shows that 



effect of the interference of the o scillator partial spectra, 
which was discussed earlier .fl^'^ 

We will consider as an example the moments (a") 
for telegraph noise. This noise takes on two values £,k 
(k = 1,2), between which it switches at random at rates 
Wi2 and W21- Respectively, in the stationary regime, 
A(n,^) = A{n,^,t) has two components, A{n,^i) and 
^('^1^2)5 which can be considered as components of a 
vector A(n); operator W becomes a 2 x 2 matrix, and 
Eq. ( 24 ) can be written as 



M(n)A(n) = inF' A{n - 1), (25) 
M(n) = n(A* + + in^c^, + W, 

where / and (T^ are the unit matrix and the Pauli matrix, 
respectively, and 



W 



W12 
-W12 



~W2l 
W21 



C = (ei+6)/2; 



= ('^i — ^2)72 is the amplitude of the frequency noise. 
Equation ([25| has a simple solution: 



(26) 



A(n) = {iF'Y n M-\k)A{0), 



functions A(n,^,t) with the same n but different ^ are 
coupled by the operator W . This is a direct analog of the 



where A(0) = {W2i/W,Wi2lW), with W ^ W12 + W21 
being the total switching rate. 

For the noise amplitude |^c| ^ W,T the absorp- 
tion spectrum of the oscillator, which is given by 
Im [A(l,^i) -(- ^(1,^2)], has two distinct peaks. On the 
other hand, for j^d ^ max F, 14^ the peaks are n ot re- 
solved and the spectrum is a single-peak curve ,^^21211 which 
makes it complicated to identify the presence of the fre- 
quency noise. The moments of the complex amplitude 
are advantageous in this respect, as discussed in the next 
section. 



VIII. COMPARING DIFFERENT TYPES OF 
FREQUENCY NOISE 

We now compare the effect of different types of fre- 
quency noise on the moments of the complex amplitude. 
We consider three common types of noise, broadband 
Gaussian and Poisson noises and telegraph noise; in the 
latter case, we choose symmetric noise with W12 = W21 
and ^2 = —£,1- In Fig. [TJ we show the dependence of the 
second moment of the complex amplitude scaled by the 
squared mean amplitude on the frequency of the driving 
field. It is obtained from Eqs. ( 12 ) and ( 26 1 . If there is no 
frequency noise, we have — 1. In the presence 

of noise the ratio | {u^) / {u)'^\ can be smaller or larger than 
one. As seen from Fig. [T] this ratio displays resonant de- 
pendence on the field frequency. It most strongly differs 
from one near resonance, where ujp = wq. As expected, 
the difference increases with the noise intensity. 
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FIG. 1: (Color online) The scaled second moment of the com- 
plex amplitude of forced vibrations u as a function of the fre- 
quency of the driving field, &lo — ujf — ^o- The solid, short- 
dash and long-dash curves show the results for the Gaus- 
sian, Poisson, and telegraph noises, respectively. The left and 
right panels refer to the relative Gaussian noise intensities 
D /T = 0.1 and 1; note the different scales. For the Poisson 
noise on these panels, we took g = 1 and the same intensities 
as the Gaussian noise, vg^ ji — D. For the telegraph noise, 
we chose W12 = W21 = D and the variance = 2Z)^. 



For weak noise, | (u^) / (u) 
tensity, and — 1 ' 



I is linear in the noise in- 
Re - is 

given by Eq. pTI). Out of the three types of the noise 



discussed in this section, only for Poisson noise, we ob- 
serve that > 1 in a certain frequency range. 
An interesting feature of this noise is that, even though 
— 0, that is, Wq is the average oscillator frequency, 
the frequency of the driving field where the field absorp- 
tion is maximal , i.e., where Im u* is maximal as a func- 
tion of oji?, is located for ujp — ujq = —v{g — sing), see 
Eq. ( [lO| ). This shift of the maximum of the absorption 
spectrum is seen in the inset in Fig. [2] 

Figures [2] and |3] show higher-order moments of the 
complex amplitude for the two values of the frequency 
detuning, Slj/T = —1 and Slj ^ (exact resonance). 
The absorption spectra of the oscillator for the chosen 
noise parameters are shown in the inset of Fig. [2] . For 
Gaussian and Poisson frequency noises, these spectra are 
Lorentzian, and therefore the presence of the noise can- 
not be inferred from the spectrum. For a telegraph noise, 
the spectrum is non-Lorentzian, but clearly is close to a 
Lorentzian curve, even though the width of the spectral 
peak is increased by a factor ^ 2. At the same time, 
the moments of the complex amplitude unambiguously 
demonstrate the presence of frequency noise. 

The dependence of the normalized moments on the or- 
der of the moment is nonmonotonic. It is very specific 
and markedly different for different types of noise. This 
is seen both in the real and imaginary parts of the mo- 
ments and in their absolute values. We do not show the 
imaginary parts of the moments at exact resonance, since 
they are small there. It is seen also that, for Gaussian 
noise, the moments decrease more rapidly than for other 
noises we study, as it would be expected from Eq. (12 1. 
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FIG. 2: (Color online) The real (upper panel) and imaginary 
(lower panel) parts of the normalized moments of the com- 
plex amplitude of forced vibrations. The squares, circles, and 
triangles show the results for the Gaussian, Poisson, and tele- 
graph noises, respectively; the lines are guides for the eye. 
The data refer to Sto/V = —1 and the same noise parameters 
as in the right panel of Fig.[l] Inset: dependence of the scaled 
complex amplitude Im {u*p} = (Aljf/F) Im(u*) on the driv- 
ing field frequency for the different types of frequency noise; 
the plotted quantity gives the oscillator absorption spectrum. 
The coding of the curves and the noise parameters are the 
same as in the main figure. 




We note that the results for the moments of the com- 
plex amplitude of classical vibrations immediately apply 



FIG. 3: (Color online) The real part and the absolute value of 
the normalized moments of the complex amplitude of forced 
vibrations at resonance, Sui = 0. The coding of the curves and 
the values of the noise parameters are the same as in Fig. [2] 



to quantum vibrations, (-u")* — (/i/2wo)"^^(a")- 

Measurements of the moments and correlators of the 
complex amplitude can be done by standard homodyne 
detection in which the in-phase and quadrature compo- 
nents of the oscillator displacement are recorded as func- 
tions of time. This procedure is standard for classical 
oscillators. For quantum oscillators, measuring the mo- 
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ments (a") is simplified by the fact tliat operators a"(t) 
with different n but the same time t commute with each 
other. The moments can be immediately found, for ex- 
ample, from the Wigner distribution, which can be mea- 
sured by means of Wigner tomography.'^- Experimental 
observation of the moments for microwave photons using 
a different procedure was reported recently.!^ 

IX. CONCLUSIONS 

The results of this paper show that, for vibrational 
modes, the presence of frequency noise can be revealed 
and the statistics of the noise can be studied using the 
moments of the complex amplitude of forced vibrations 
(«,") oc (a")*. The moments can be directly measured in 
the experiment. In the presence of frequency noise, they 
differ from the corresponding powers of the average com- 
plex amplitude (u)" oc ((a)")*. The moments display a 
characteristic dependence on the frequency of the driv- 
ing field and the moment number and are very sensitive 
to the noise statistics. This is illustrated using as exam- 



ples Gaussian and Poisson noises with bandwidth that 
significantly exceeds the oscillator decay rate, as well as 
a telegraph noise. 

Explicit expressions are obtained for the moments of 
the complex amplitude in the case of broadband noise 
with arbitrary statistics. A general formulation is de- 
veloped for Markov noise, which reduces the problem of 
calculating the moments to a set of linear equations. Ex- 
plicit results for the variance and the third cumulant of 
the complex amplitude are obtained also for an arbitrary 
noise provided the noise is weak; the third cumulant of 
the amplitude in this case is proportional to the third cu- 
mulant of the noise. Even for weak frequency noise, the 
proposed method allows revealing it irrespective of the 
intensity of additive noise in the oscillator. The results 
apply to both classical and quantum oscillators. 
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